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1 Abstract. We investigate the accretion onto luminous bodies with hard surfaces within the 

CNj ' framework of newtonian theory. The accreting gases are assumed to be polytropes and their 

selfgravitation is included. A remarkable feature of the model is that under proper boundary 
conditions some parameters of the sonic point are the same as in the Bondi model and the 
relation between luminosity and the gas abundance reduces to an algebraic relation. All that 
holds under assumptions of stationarity and spherical symmetric. Assuming data that are 
required for the complete specification of the system, one finds that generically for a given 
luminosity there exist two solutions with different compact cores. 
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1. Introduction 

Stationary accretion of spherically symmetric fluids with luminosities close to the Eddington 
limit has been investigated since the pioneering work of Shakura pQ. In the gas pressure and 
its selfgravity have been ignored and the analysis has been purely newtonian. Later researchers 
included both these aspects [2] and extended the analysis onto relativistic systems 

This paper deals with the complete newtonian picture (which we will refer to as the Shakura 
model), in which pressure and selfgravity are included. It is assumed that a ball of a polytropic 
selfgravitating gas accretes stationarily onto a body with a hard surface; both are spherically 
symmetric. Our goal is to study a kind of an inverse problem. Let be given the luminosity, 
total mass, asymptotic temperature and the equation of state of the gas. Assume that the 
gravitational potential has a fixed value at the boundary of a compact body but the radius of the 
body is arbitrary to some extent. One can easily check that these data are sufficient to determine 
solutions of the Shakura model; any additional information would lead to a contradiction. On 
the other hand, it is well known that nonlinear problems can exhibit nonuniqueness. In our 
concrete example the nonuniqueness would manifest itself in multiple solutions which have 
different central compact cores. That issue is interesting in itself within the framework of 
astrophysics. The fundamental question is whether observations can distinguish compact bodies 
endowed with a hard surface (say gravastars [3]) from black holes. Narayan presented 
observational arguments in favour of the existence of black holes. Abramowicz, Kluzniak and 
Lasota [Sj raised several objections, pointing out that present accretion disk models are not 
capable to make distinction between accretion onto a black hole or a gravastar. 

The Shakura model is the simplest selfcontained system that can be interpreted as a radiating 
system. Since we work at the newtonian level, the issue of making distinction between a black 



hole and a compact body with a hard surface is not addressed. Our results show, however, that 
there is an ambiguity. There can exists at least two systems (two compact stars with a hard 
surface, having different masses) to the inverse problem formulated above. This is consistent with 
the recent general-relativistic analysis IB] , valid in the newtonian limit (with a reservation) , 
which shows that the mass accretion rate M behaves like y 2 (l — y). Here y = M*/M is the ratio 
of (roughly) the central mass to the total mass. Therefore the maximum of M is achieved at 
y = 2/3. The luminosity is equal to the mass accretion rate M multiplied by the available energy 
per unit mass, the potential 4>(Ro), at the hard surface of a compact star. Since the boundary 
value <ft(Ro) is fixed, this would mean that there appear two weakly luminous regimes: one rich 
in fluid with rrif ~ M and the other with a small amount of fluid, rrif/M <C 1. That would 
suggest a nonuniqueness — given a luminosity x, one would have two systems with different 
y's. The situation in the Shakura model is not so simple, because the luminosity impacts also 
the accretion rate M and there emerges a complex functional relation y = y(x). Our analysis 
shows that under suitable boundary conditions this relation y = y(x) is well approximated by 
an algebraic equation, and then one shows that the nonuniqueness again emerges in the Shakura 
model. 



2. Main equations 

We will study a spherically symmetric stationary flow of a fluid onto a spherical compact body. 
The areal velocity U of a comoving particle labeled by coordinates (r, t) is given by U (r, t) = dtR, 
where t is a comoving (Lagrangian) time, p denotes pressure, L(R) and Le are the local and 
Eddington luminosities, respectively. The mass m{R) is given by d R m{R) = 4nR 2 g. The mass 
accretion rate reads 

M = -4ttR 2 qU, (1) 

where g is the baryonic mass density. Since the fluid is a polytrope, we have the equation 
of state p = Kg r with a constant T, (1 < T < 5/3). The steadiness of the collapse means 
that all characteristics of the fluid are constant at a fixed areal radius R: dtX\ji =cons t = 
(dt — (dtR)dji)X = 0, where X = g,U,a 2 . The speed of sound is given by a = ^J'd s p. Strictly 
saying, a stationary accretion must lead to to the increase of the central mass. This in turn 
means that the notion of the "steady accretion" is approximate — it demands that the mass 
accretion rate is small and the time scale is short enough, so that the quasilocal mass m(R) does 
not change significantly. The total mass m(R 00 ) will be denoted by M. It is assumed that the 
radius Roq of the ball of fluid and other boundary data are such that 

^«^M <<<4 . (2) 

-too 

The value of the gravitational potential at the surface of the compact body is fixed to be —<j)o- 
The steadily accreting polytropic fluid is described by a system of nonlinear integro- 
differential equations. They consist of the Euler equation 

Gm(R) l n L(r) 

R = — ^r 1 - -°rp + a -^> ( 3 ) 

the mass conservation 

d R M = 0, (4) 

and the energy conservation 

■fa 2 a 2 U 2 \ 

L - L(r) =M j^-^-^- m ■ (5) 



The constant a is given by a = ut I ^vn v c, Lq is the total luminosity and cft(R) is the newtonian 
gravitational potential. We neglected in Eq. (J5J) the term — jr- since it is much smaller than 

jrfHr, due to the boundary conditions that have been displayed above. 

The Eddington luminosity is calculated here for the whole system, that is Le = GM/a, while 
the total luminosity Lq is equal to the product of the mass accretion rate and (j)Q = \(f)(Ro)\, 
where i?o is an areal radius of the boundary of the compact body. Lq = McpQ pQ. The quantity 
0o is fixed, as said above, and Rq in the constructed configurations is by definition the radius at 
which the absolute value of the potential equals to 4>q. Rq can be found only after the solution of 
the whole system is constructed. After some algebra (which consists of differentiating Q with 
respect to R, using the equality djia 2 /(T — 1) = a 2 djilng = djip/g and combining the obtained 
equation with Q) one arrives at 

d R lnL = ^-. (6) 

This can be immediately solved, 

Here appears a new quantity Rq = GM/\(J)(Rq)\. Rq represents a kind of the size measure of the 
compact body. For test fluids (when the absolute value of the potential is just the product of 
the gravitational constant G by the mass divided by the radius) one has Rq = Rq, while in the 
general case R > Rq. 

One can prove, using the same arguments that Bondi applied to the accreting newtonian 
stars with test gases that there exists a unique solution in the case without selfgravitation. 
We believe that this is just a matter of a simple (and perhaps tedious) work to prove analytically 
the existence of solutions of the Shakura model. In any case, there exist numerical solutions 11 . 



3. The mass accretion rate 

Below we assume that the accretion flow possesses the so-called sonic point, but in our opinion 
the same conclusions should be valid for flows without sonic points. The standard terminology — 
critical flows (if there is a sonic point) and subcritical flows (if the sonic point is absent) — will 
not be used in this paper. The notions "critical solutions" and "subcritical solutions" will 
appear, but in a different meaning, defined later. 

The sonic horizon (sonic point) is at a location where \U\ = a. In the following we will 
denote by the asterisk all values referring to the sonic points, e.g., a*, U*, etc. Differentiation 
with respect to the areal radius will be denoted as prime '. The mass conservation equation 
yields 

U' = -U {q'/q + 2/R) . (8) 
Inserting Q into Eq. © one arrives at 

Using the definition of the sonic point one discovers that its three characteristics, a*, f7* and 
M*/i?* must satisfy following relations 

2 tt9 GM* A L*a \ GM* ( H L*M 



< - ' = -YiC V - ok) = w: I 1 " Y^m, ■ < 10 » 



Eq. (|10|) implies that if there exists a sonic point then 

L*M 



L E M* 

The velocity U and the mass density g can be expressed as follows 

1/(T-1) 

and 



< I- (11) 



U = VM(±\ . (121 



Q = Qoo (a/a OD ) 2/{r x) . (13) 

Here [/* is the negative square root and the constant is the asymptotic mass density of a 
collapsing fluid. Two of the four classes of stationary solutions that cross at the sonic point can 
be assigned physical meaning of an accretion or a wind. 
Let us introduce auxiliary variables 



L .._Af, _ Rq 

Tj<a 



L - L* - 2Mal 



yL E - L* 



X GM * tlA\ 

The necessary condition pip for the existence of a sonic point can be now formulated as 

xexp(— x'y) < y. (15) 
A straightforward calculation yields the following equation 

al ~ = (T - 1) - . (16) 



M 



We shall assume that 7 < 1 and j;7<1. 
Eq. (JTU|) implies that 



2 GM* ( ^ x C\ GM* ^ x 



aj = ^1 1 - - exp (-x 7 ) « 1 - - . (17) 

One obtains from p4p (keeping only the first term in the expansion of the luminosity function L*) 

GM x MGM* 



Ln — L* — xLr, 



\R*<j)(Ro)\ y R* 



x 



2Ma . (18) 
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In the same approximation the third term in the expression for A* reads — 2Ma^^j; therefore 
A*/M ~ 0. One can show, using a reasoning similar to that provided in the next section, that 
\£* = 4G-7T /^°° drrg is much smaller than a^. 

Thence the sonic point formula p6p simplifies to 



That is one of the main results of this presentation. The reader can be surprised by the fact that 
we obtained the same result as in the Bondi model, in which the luminosity and selfgravitation 
are absent ^0]- The explanation for that lies in our boundary conditions (J2J). Relaxing them 
(for instance assuming a nonzero asymptotic falloff velocity U) would invalidate formula (|19j) 
and also the remaining results of this work. We would also like to warn the reader that it is 
assumed in the above derivation that T is somewhat smaller than 5/3. 



4. Luminosity 

In this section we return to the question stated at the beginning. Namely, let be given the total 
mass M, the size Roo, the luminosity Lq and the asymptotic temperature T (or equivalently, 
the asymptotic speed of sound a 2 ^). We show that in the generic case there exist two accreting 
systems. 

The rate of the mass accretion M (see Eq. ([T])) can be formulated in terms of the 
characteristics of the sonic point. After some algebra one derives from Eqs (jlO|) and (|T|) 

(5-3r) 

M = G 2 7rM 2 ^(y-xexp(-x 1 )) 2 { ) ^ . (20) 

Now, one can can show in a way similar to that applied in the case of relativistic accretion 
|Hj that under the previously assumed conditions and ary <C 1 one has g^ — Xoo 

(M - M*) = 

Mxoc (1 — y) for T G (1, 5/3 — 5) with some small 5. 

We shall describe briefly main stages of this calculation. It is easy to see, that Eq. © yields 

a 2 (R) r-l,, n , 

— ^ < 1 - — 2-tf(i2). (21) 

Now observe that \(f)\ < GM/R; therefore 

a 2 (R) _ , GM 

_L2<i + (r-i) (22) 

Returning now to the expression of the mass density given in (|13|) and replacing the ratio 

a 2 (R)/al Q by the right hand side of <0, one arrives at g < g^ (l + (T - 1) ^t-) 1 1 • 
Applying this bound in M — M* = J v dV g and invoking the asymptotic conditions © one 
obtains the needed approximation — the equality g^ Mxco (1 — y)- The proportionality 
constant Xoo is roughly the inverse of the volume of the gas outside of the sonic sphere. 

Now notice that the total luminosity Lq = (J)(Rq)M. Replacing in M (see Eq. (|20(l ) the ratio 
a 2 1 a 2 ^ by the right hand side of ()19j) , one obtains 

o (5-3r) 
2 M 3 2 ( 2 A 

L = cj) G nxoc-^- (1 - y) (y - x exp (-xj)) ^ - _ 3r J . (23) 

It is convenient to rewrite this equation using the relative luminosity x = Lq/Le- 

x = (3{l-y){y-xex.p{-x"f)) 2 . (24) 



Here (3 is a dimensionless numerical factor depending only on the boundary data characterizing 

2 , . (5-3r) 

the system in question, (3 = Xoo&ootGTr ( r 2 v ) 2(r_1) . We have assumed that the luminosity 



x is given and our aim is to find all possible solutions y(x). It appears useful to start from a 
critical point and then show the existence of a bifurcation. Let us define 



F(x, y) = x - (3 (1 - y) (y - x exp (-X7)) 2 . (25) 

Let (a, 6) be a zero of F, that is F(a, b) = 0. We will say that (a, b) is critical if d y F(x, y)\ a ,b = 0. 

One can show that there exist at least two solutions y(x, (5) for any parameter (3 (0 < 7 < 1) 
and for the relative luminosity x smaller than a. The obtained results are following. 

Theorem 

(i) There exists a unique critical point x = a, y = b of F, with a and b satisfying bounds 
< a < b < 1 and the relation b = 2+aex | ( ~ a7) . 

(ii) For any < x < a there exist two solutions y(x)l bifurcating from (a, b). They are locally 
approximated by formulae 



y + = b± V( a - x)(b + aexp(-aj)(l - 2a 7 )) 
\/f3(b — aexp(— «7))(1 — aexp(— 07)) 

(iii) The relative luminosity x is extremized at the critical point (a,b). 
Proof 

The two conditions F = and d y F = yield two equations 

a - - 6) (6 - aexp(-a7)) 2 = 0, 

£>-aexp(-a 7 ) =2(1-6). (27) 

From Eqs (|27|) one immediately obtains 

b= 2 + aexp(-a7)_ (2§) 
3 

Let us remark that from the above one can get also 

a = 4/3(1 - 6) 3 ; (29) 
inserting that into the second equation in (|27|). one arrives at 

6 = I + 1r (1 " 6)3 exp ( _4/?(1 " 6)37 ) • (30) 

Both sides of this equation are continuous functions of b and at b = the right hand side of (|30j) 
is bigger than the left hand side, while at b = 1 the opposite holds true. That guarantees the 
existence of a solution. Closer investigation allows one to conclude that the critical solution is 
unique. This fact can be used in order to show that solutions y(x) bifurcating from (a, b) extend 
onto the whole interval x G (0,a). The argument relies on the use of the implicit function 
theorem and the uniqueness of the critical point (a,b). The specific form of a solution close to 
a critical point can be obtained in a standard way. Put x = a + e, y = b + y(e) into F(x, y) = 
and expand F keeping a few terms of lowest order. One obtains a reduced algebraic equation 
(known in the mathematical literature as the Lyapunov- Schmidt equation) 

(1 + 2(3 (1 — b) {b — aexp(— 07)) (aexp(— 07)7 — exp(— 07))) ex 

P (36 - 1 - 2a exp(-a 7 )) y{ef = 0. (31) 



At the critical point 36 — 1 — 2aexp(— 07) = 1 — aexp(— 07) (see the second of Eqs (|27jl) while 



2(3 (1 — b) (b — aexp(— 07)) (aexp(— 07)7 — exp(— 07)) = 

= 2a (aexp(— 07)7 — exp(— 07)) /(b — aexp(— 07)). (32) 

Inserting (|32[) into ()31|) and finding y(e) from the latter immediately leads to the approximate 
solution of (|n|). 

The proof of the third part of the Theorem goes as follows. Let (xo,yo) be a non-critical 
solution of F(x, y) = with the domain belonging to a subset xe~ x < y of the square < x < 1, 
< y < 1. The implicit function theorem guarantees the existence of a curve x{y) such that 
F(x(y),y) = for y belonging to some vicinity of yn. Along this curve one has 

dx = (y- xe- x ~f) (% - 2 - are"* 7 ) , . 

dy 1 + 2/3(1 - y)(7/ - xe" :r T)(e- a; T -7xe- x T)' 1 ' 

Notice that the denominator is strictly positive while the nominator vanishes only at critical 
points. That proves the assertion. Now it is clear from the form of approximate solutions 
constructed in part (JnJ) that bifurcating solutions exist only if x < a (our bifurcation is called 
sub critical in the mathematical literature); therefore x = a is the extremal value of the relative 
luminosity. That ends the proof. 

The bifurcating solutions y(x) are shown in the Figure, for the parameter (3 = 1000 and 
under the assumption 7 = 0. Notice that coordinates (a, b) of the critical point increase 
with the increase of the parameter (3. In order to show that, observe that Eq. ()28|) implies 

Bifurcation point = 0.820, 0.9236003642 
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Figure 1. The mass abundance is put on the ordinate while the relative luminosity x is shown 
on the abscissa. The circle encloses the bifurcation point (a,b). 



4f = e a7 (l — 07). This expression is bigger than zero, since 07 < 1. Differentiating Eq. (|29j) 
with respect to (3 gives 4| = 4(1 - b) 3 /(^ + 12/?(1 - 6) 2 ) > 0; this inequality is due to the 
positivity of Now we have along the critical curve % = gf therefore also t§ > 0. Let 
us remark that if the parameter 7 vanishes then one can explicitly find b by solving ()3U|) , The 
solution reads 

(/3 2 + /? 3 / 2 v / r+?) 2/3 -/? 
^ = 1 tr ; (34) 

(129)) gives then a(/3) and one can check explicitly that both a and & monotonically increase with 
the increase of f3. 

It is of interest that at critical points the parameter b is not smaller than 2/3. From the 
analysis of the foregoing formulae it is clear that b is close to the value 2/3 if (3 is small, i.e., 
when the relative luminosity a is small (notice that a < (3). The value 2/3 appears in the study 
of the mass accretion in irradiating systems 0. The fluid abundance is equal to 1 — b; thus 
we can infer that critical configurations have less fluid than 1/3 of the total mass, and that the 
upper bound 1/3 is achieved in the limit of vanishing radiation. The maximal mass accretion 
rate does not correspond to y = 2/3, as in systems with no radiation 0, but to a somewhat 
larger value. 



5. Conclusions 

In this presentation we have assumed the existence of an accreting system satisfying particular 
boundary conditions. We would like to stress out that (jij) the detailed behaviour of the accreting 
flow is given by the integro-differential nonlinear Eqs ©-© and (jnj) due to our boundary 
conditions some essential features can be deduced from the algebraic equation H24|) . The question 
arises whether this projection of the original equations onto the algebraic problem in fact takes 
place. We checked numerically that there do exist appropriate solutions and that the relative 
error made in the above approximations is of the order of 10 -3 . An explicit example can be 
found in [TT] . 

It is illuminating to repeat the preceding discussion in the case when 7 < 1. A number of 
simplification occurs now, since one can approximate the term xexp(— ary) x. From Eqs (|29j) 
and (|34j) one can see that a can be as close as one wishes to 1 for sufficiently (3. Expressing that 
conclusion in more intuitive terms: the total luminosity of the accreting system Lq approaches 
the Eddington luminosity Le if the parameter (3 goes to infinity. In particular, from the 
expression for (3, if the product (M/Mo) 2 /(aoo/aooo) 3 (where Mo and aooo are some reference 
quantities) is large enough, then Lq is close to Le- At the critical point there is only one 
solution. The two bifurcation solutions are characterized by the luminosity x < a and their 
(central) mass parameters y± and 2/2 can be markedly different only for x <C a. This can be seen 
from the point (JmJ) of the Theorem and confirmed by the Figure, which shows the branching of 
the two solutions from a critical point in a class of examples. If a given system radiates with 
the luminosity close to its critical luminosity (and, in particular, to the Eddington limit), then 
the cores corresponding to the bifurcating solutions have similar masses. The necessary and 
sufficient conditions for having two accreting compact stars with hard surfaces that have masses 
satisfying condition Mi <C M2 is that the luminosity is much smaller than Le or the critical 
luminosity o, respectively. 
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